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SECTION  1 
INTRODUCTION 


This  report  presents  scaled  numerical  solutions  for  the  one  dim- 
ensional time-dependent  SGEMP  electron  emission  boundary  layer  problem  In  the 
two  cases  in  which  the  incident  X-ray  flux  is  either  linearly  rising  in  time 
(a  linear  ramp)  or  is  constant  in  time  starting  at  t>0  (a  step  function) . 
Two  electron  energy  spectra  are  considered  in  each  case,  an  exponential 
spectrum,  and  a linear-times-exponential  spectrum.  The  electrons  are  as- 
sumed to  be  emitted  with  a cos6  angular  distribution.  The  scaled  solu- 
tions given  here  provide  the  solutions  to  real  problems  with  these  energy 
spectra  for  arbitrary  X-ray  flux  rise  rates  (for  the  linear  ramp)  or  arbi- 
trary fluxes  (for  the  step  function) . 

Most  X-ray  pulses  start  off  rising  linearly  in  time.  In  many 
cases,  especially  at  high  fluences  or  long  rise  times,  much  of  the  interest- 
ing boundary  layer  dynamics  is  taking  place  while  the  pulse  is  still  linearly 
rising.  Since  a single  solution  suffices  for  all  rise  rates  during  this 
portion  of  the  pulse,  it  is  worthwhile  to  present  this  solution. 

In  addition,  if  tho  pulse  rise  time  is  very  short  compared  to  an 
electron  turnarv-und  .v'’''  and  compared  to  the  subsequent  time  over  which  the 
remainder  of  the  pulse  var;)  s significantly,  then  the  bounder'-  ’ ‘vn.l 
behave  approxx  tely  as  if  the  ^-ray  pulse  were  a step  function.  We  also 
present  the  solu  on  for  this  case. 


The  scaled  e...  ^'ous  controlling  the  boundary  layer  were  pi.,  -nted 
in  Reference  1.  There  it  was  snow,  that  if  x-ray  pulse  rises  in  time 


' Longmire,  C.  L. , and  N.  J.  Carron,  Scaling  of  the  Time- Dependent  SGEMP 
Boundary  Layer,  Mission  Research  Corporation,  MRC-R-262,  DNA  3975T, 
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April  1976  (U) 


like  a power  of  time 


F - Rp  tP  (1) 

where 

0 i p < « , (2) 

F is  the  incident  flux  (cal/cm  /sec),  and  Rp  a constant,  then  the  solu- 
tion scales  completely.  For  a given  electron  energy  spectrum  and  given  p, 
a single  solution  suffices  for  all  Rp.  The  solutions  given  here  correspond 
to  p ■ 1 (Rj^  is  the  flux  rise  rate)  and  P“0  (R^  is  the  constant  flux). 

The  solution  is  a one- dimensional  one  and  should  be  meaningful 
out  to  distances  small  compared  to  the  lateral  dimensions  of  the  target. 

We  give  plots  of  electric  field,  potential,  number  density,  and 
current  density  as  a function  of  position  and  time.  These  are  presented 
both  as  a function  of  x (distance  from  the  target  surface)  for  various 
times,  and  as  a function  of  time  for  various  x.  in  addition  the  dipole 
moment  per  unit  area  of  the  layer,  and  its  first  and  second  time  derivatives 
are  given  as  a function  of  time.  The  spatial  integral  of  the  current  density 
is  the  same  as  the  time  derivative  of  the  dipole  moment. 

The  current  densities  given  here  may  be  useful  as  prescribed  cur- 
rents in  Maxwell  solving  codes.  The  second  time  derivative  of  the  dipole 
moment  may  be  useful  for  estimating  dipole  radiation.  The  dipole  moment 
may  be  useful  for  estimating  distant  quasi-static  dipole  fields. 
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SECTION  2 

COMPUTATION  METHOD 


The  numerical  solutions  were  obtained  with  the  one  dimensional 
particle  moving  code  SCALID.^  Time  and  length  scales  are  determined  by 
the  naturally  occurring  electron  plasma  properties.  Times  are  measured 
as  multiples  of  an  effective  plasma  period  T^  given  by 


where  m(e)  is  the  electron  mass  (charge)  and  is  an  average  surface 
number  density,  and  lengths  are  measured  as  multiples  of  an  effective  plasma 
Debye  length 

A - V Tp  (4) 

where  v is  the  average  normal  velf^  ity  of  the  emitted  electrons.  These 
times  and  lengths  are  scaled  out  before  beginning  computation,  so  computed 
times  and  lengths  are  on  the  order  of  unity.  Similarly  the  dependent  dy- 
namical variables  have  a natural  magnitude  which  is  factored  out. 

If  the  pulse  is  linearly  rising,  Nj  in  Equation  (3)  is  itself 
time  dependent,  so  that  Equation  (3)  as  it  stands  does  not  define  a time 
unit.  In  the  general  case.  Reference  1 shows  that  a plasma  period  can  be 
taken  to  be 

* Carron,  N.  J. , Description  of  the  Code  SCALID  for  Calculating  the  One- 
Dimensional  SGEMP  Boundary  Layer.  Mission  Research  Corporation,  MRC-R-267, 
May  1976  (U) . 
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(5) 


m V 


47re‘Y  R 


1 

2+p 


where  Y is  the  material  yield  (elec/cal) 


We  assume  the  emitted  electrons  have  a cos6  angular  distribu- 
tion, where  0 is  the  polar  angle  measured  from  the  normal.  The  shapes 
of  the  two  energy  spectra  we  consider  are  shown  in  Figure  1.  If  the  elec 
tron  energy  spectrum  is  exponential 


dn  1 E/E 


(elec/keV) 


C6) 


where  E is  the  electron  energy  and  E is  the  exponentiation  energy,  then 
the  average  normal  velocity  is 

If  the  spectrum  is  linear-times-exponential 


dn  ^ E „-E/E  , , 

■jn"  ® * (elec/keV) 

dE  -^2 

then 


(8) 


v 


(9) 


In  the  code  we  have  chosen  a very  small  time  step,  i 0.02  T . At 
all  interesting  times  there  were  several  thousand  particles  in  the  first 
few  Debye  lengths.  The  particles  were  given  a finite  length  of  0.2  Debye 
length  for  the  purpose  of  computing  the  charge  density  and  current  density. 
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Energy  Spectrum 


Figure  1.  Exponential,  and  llnear-times-exponentlal  energy  spectra. 
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These  quantities  were  further  averaged  over  three  spatial  cells  with  a sim- 
ple weighting  scheme.  The  coll  size  was  -0.1  Debye  length,  and  there  were 
always  several  hundred  particles  in  each  cell. 

The  current  density  J(x,t)  was  computed  as  the  sum  of  particle 
weights  times  their  velocities,  rather  than  as  3E/3t  where  E is  the 
electric  field.  At  later  times  when  J is  small  and  due  to  approximately 
equal  numbers  of  particles  with  positive  and  negative  velocities,  irrtigu- 
larities  and  oscillations  occurred  whose  precise  value  was  dependent  on 
the  time  step,  emission  particle  weighting  scheme,  etc.  In  the  plots  pre- 
sented here  we  show  an  average  smooth  dashed  line  through  these  unreliable 
values.  The  plots  of  current  density  vs.  time  at  various  x were  some- 
times constructed  from  3E/3t  when  this  method  was  more  reliable. 
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SECTION  3 

LINEARLY  RISING  PULSE.  EXPONENTIAL  SPECTRUM 


In  this  section  we  present  results  for  a linearly  rising  time 
history  Cp"l)  and  for  the  exponential  energy  spectrum,  Equation  (6). 

In  the  following  equations  we  give  dimensional  units  for  the 
dynamical  variables.  When  numbers  are  given  in  the  equations,  we  use  these 
units  for  E,  Y,  and  : 

Electron  exponentiation  energy  1"  (keV) 

Yield  Y (elec/cal) 

2 2 

X-ray  rise  rate  (cal/cm  /sec  ) 


Equations  (5)  and  (7)  show  the  time  unit  to  be 


and  the  unit  of  length,  Equation  (4),  is 


CIO) 


A - 7.796  X 10 


^2  11/3 


y Ri 

^ j 


cm 


Cll) 


The  units  of  electric  field  (E^)  and  potential  (<I>^)  are  taken  to  be 
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'1  e T 


- 8.955  X 10 


* ^ C 2^5 

*1  ■ ■ T ? 


-5  Ty 


[e  y r,]- 


Voits/m  , 


698  E Volts 


The  units  of  number  density  (Np  and  current  density 

VP  rp 

* ■'1  in  * ‘'i  \ 

N,  - - 6.349  X ^ cm  , (14) 

^ V E 

Jj  - e Nj  V • 1.126  X 10"^^ [Vi"  Amps/in^  . (15) 

The  units  of  dipolo  moment  per  unit  area  (P^)  and  its  time  derivative; 
(Pp  are 

‘’l  “ ^ ^ • 187 


e 6.181  X 10  E CoulM 


P^  » ^ - 8.786  X 10’®  [e^^^  Y Amps/m 


Finally,  the  second  time  derivative  of  the  dipole  moment  per  unit  area  has 


p’^  = -5“  = 1.249  X io“®  j^E  Y Amps/m/sec 
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For  example,  if  E ■ 5 keV,  the  material  yield  is  Y ■ 5 x 10 
elec/cal,  and  the  flux  is  rising  at  a rate  R, 


12 


10^^  cal/cm^/sec^,  then 


T ■ 2.50  ns 
P 

X ■ 6.19  cm 

« 5.64  X 10^  V/m 
- 3.49  kV 
■ 5.04  X 10^  cm“^ 

Jj  ■ 2.00  X 10^  Amps/m^ 

Pj  - 3.09  X 10"®  Coul/m 

Pj  “ 12.38  Amps/m 

• # 9 . 

Pj  B 4.96  X 10  Amps/m/sec 


Figures  2 through  9 show  the  electric  field,  potential,  number 
density,  and  current  density  as  a function  of  distance  from  the  surface  and 
time,  scaled  to  X and  and  Figures  10  through  12  show  the  layer's 

dipole  moment  per  unit  area  and  its  first  two  time  derivatives.  The  units 
are  those  of  Equations  (10)  through  (18) . 


11 


2 


I 


1 


f 


0 

0 1 2 3 4 5 

x/X 


6 

5 
4 

LU 

uj  3 

2 
1 
0 

0 2 4 6 8 10  12  14  16  18  20 

x/X 

Figure  2.  Electric  field  vs.  x at  various  times. 
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Figure  3.  Potential  vs.  x at  various  times. 
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SECTION  4 

LINEARLY  RISING  PULSE,  LINEAR-TIMES- 
EXPONENTIAL  ENERGY  SPECTRUM 


This  section  presents  results  for  a linearly  rising  time  history 

(p>l)  and  for  the  llnear-tlmes-exponentlal  energy  spectrum,  Equation  (8) . 

¥,  Y,  and  R,  are  In  the  same  units  as  In  the  previous  section  (keV,  cal~\ 
2 2 

cal/cm  /sec  ] . 


Equations  CS)  and  (9)  show  the  time  unit  to  be 


r .Lj ViJ 

P L4^e2  YrJ 


Tll/3 


- 0.8053 

and  the  unit  of  length  Is 


¥*11/3 


sec 


■ V T - [ 

P Lie  e 


^2-11/3 


2 Y R, 
m ] 


1.339  X 10‘ 


cm 


The  units  of  electric  field  (E^  and  potential  ($j^)  are 


_ m _v 
“ e T. 


[2tt^  ^ E Y 

tr 

1.173  * 10'^  Y rJ^^^  Volts/i 


(19) 


(20) 


(21) 


20 


(22) 


1 


rr  f 

2 e 


1.571  X 10^  E 


Volts 


The  units  of  number  density  (N^)  and  current  density  (J^)  are 


N 


1 


v 


The  units  of  dipole 
(fj)  are 


4.845  X io“^° 
e Nj  V 

1.290  X 10"^®  [Vf”  Amps/ui^  (24) 

moment  per  unit  area  (P^)  and  its  time  derivative 


B 


1/3 


cm 


-3 


(23) 


fi  ■ • "i  ■ A 

■ 1.391  X 10"®  E Coul/m  , (25) 

^ . 1.727  X io‘®  [e®^^  Y Amps/m  (26) 

P 

The  units  of  second  time  derivative  of  the  dipole  moment  per  unit  area  is 

p’j  ■ ■ 2.144  X 10"®  [e  Y rJ^^^  Amps/m/sec  . (27) 

P 

Figures  13  through  23  show  the  dynamical  solutions  for  the  case  of  a linearly 
rising  time  history  and  a llnear-times-exponentlal  energy  spectrum,  using  the 
units  of  Equations  (19)  through  (27). 
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Figure  13.  t:]ectr1c  field  vs.  x at  various  times 


SECTION  S 

CONSTANT  PULSE,  EXPONENTIAL  SPECTRUM 


In  this  section  we  present  results  for  a constant  X-ray  pulse 
(p«0)  and  for  the  exponential  energy  spectrum,  Equation  (6) . 


We  give  here  the  dimensional  units  for  the  dynamical  variables. 
When  numbers  are  given  we  use  these  units  for  Y,  and  R^: 


Electron  exponentiation  energy 
Yield 

X-ray  flux 


E (keV) 

Y (elec/cal) 

2 

(cal /cm  /sec)  . 


Equations 


(5)  and  (7)  show  the  time  unit  to  be 


and  the  unit  of  length  Is 


(28) 


(29) 


30 


The  units  of  electric  field  (B^)  and  potential  (4^^)  are 


Ej  ■ 2-^  . 1.068  X lo'^iyfy  V/m  , (30) 

*1  ■ X El  ■ - 698  E Volts  . (31) 

The  ur.its  of  number  density  (N^)  and  current  density  (J^)  are 

N,  ■ ■ 9.025  X io"^°  cm‘^  , (32) 

^ V yjE 

■ e V - 1.602  X lo"^®  Y Amps/m^  (33) 


The  units  of  dipole  moment  per  unit  area  (P^  and  its  time  derivative 
(Pj)  are 


'i  ■ • ■ ih 


6.181  X lo‘®  E Coul/m 


1 T. 


- ■ 1.047  X 10"®  Y Rq]^^^  Amps/ 


m 


(34) 

(35) 


The  second  time  derivative  of  the  dipole  moment  per  unit  area  (Pp  has 
unit 


V,  m • 1.775  X 10"®VT 


1 T_ 


Y Rp  Amps/m/sec 


(36) 
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— 13 

Eor  example  If  B > 5 keV,  Y ■ 10  elec/cal,  and  the  flux  la 
S 2 

■ 10  cal/cm  /aec,  then 


0.88 

ns 

X 

2.19 

em 

1.60 

6 

> 

U) 

o 

X 

♦l 

3.49 

kV 

4.04 

X 10®  cm"® 

•'l 

1.60 

3 2 

X 10  Amps/m 

'l 

3.09 

X 10"®  Coul/m 

■ 35.01  Amps/m 

■ 3.97 

X 10^®  Amps/m/sec  . 

Figures  24  through  34  give  the  dynamical  solutions  'or  the  case 
of  a constant  X-ray  pulse  starting  at  t«0  and  an  exponential  electron 
energy  spectrum.  The  units  are  those  of  Equations  (28)  through  (36). 
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Figure  28.  Electric  field  vs.  time  at  various  x. 


SECTION  6 

CONSTANT  PULSE*  LINEAR-TIMES-EXPOMENTIAL 
ENERGY  SPECTRUM 


Here  we  present  results  for  the  constant  pulse  (p*0)  and  linear- 

tlmes-exponentlal  energy  spectrum*  Equation  (8).  E*  Y,  and  R are  in  the 

-1  2 ® 

same  units  as  the  previous  section  (keV,  cal  * 6al/cm  /sec). 


Equations  (S)  and  (9)  show  the  time  unit  is 

ll/2 


m V 


4Tre‘‘  Y R. 


■ 0.7227 

and  the  unit  of  length  is 


1/2 


sec 


X - V Tp  - 1.201  X 10 


K] 


1/2 


cm 


The  units  of  electric  field  (Ej)  and  potential  (♦, ) are 


'1  e T 


^ - 1.308  X 


'1- 


[VY  V R J 


1/2 


V/ra 


(37) 


(38) 


(39) 


\ c . Tf  E 


1.571  X 10^  E Volts 


(40) 


41 


Hie  units  of  number  density  (N^^)  and  current  density  (J^)  are 

V R. 


N 


1 


° 6.016  X 10‘^° 


^-15 


w 

cm"®  , 

(41) 

2 

Amps/m 

(42) 

The  units  of  dipole  moment  por  unit  area  (P.)  and  its  time  derivative 

• ^ 


(Pj)  are 


- 1.391  X io“®  B Coul/tt 


(43) 


Pj  ■ ^ ■ 1.924  X io“®  Y Amps/m 


(44) 


Finally,  the  unit  of  the  second  time  derivative  of  the  dipole  moment  per 
unit  area  is 


Vj  • ^ • 2.663  X 10"® y Pq  Amps/m/sec 
P 


(45) 


Pigxires  35  through  45  give  the  dynamical  solutions  for  the  case 
of  a constant  X-ray  pulse  starting  at  t^O  and  a linear-tlmes-exponential 
energy  spectrum.  The  units  are  those  of  Equations  (37)  through  (45) . 
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T 


Eq.  39 
Eq.  38  - 

Eq.  37 
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density  vs.  time  at  various  x. 
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ATTN:  DRXDO'RCC,  John  A.  Rosado 
AHN:  DRXDO-NP 

Commander 

Pleatinny  Arsenal 

ATTN:  SMUPA 
ATTN:  SARPA 

Comnandar 

Redstone  Scientific  Information  Ctr. 

U.S.  Amy  Missile  Command 

ATTN:  Chief,  Documents 


Chief 

U.S.  Amy  Comunlcatlons  Syi.  Agency 
ATTN:  SCCM<AD-SV.  Library 

Commander 

U.S.  Amy  Electronics  Command 
ATTN:  ORSEL 

Commander 

U.S.  Amy  Foreign  Sclanca  & Tech.  Ctr. 
ATTN:  ORXST-ISI 

DEPARTMENT  OF  THE  NAVY 

Chief  of  Naval  Operations 
ATTN:  Code  604C3 

Chief  of  Naval  Research 

ATTN:  Henry  Mullaney,  Code  427 


Director 

Naval  Reiearch  Laboratory 
ATTN:  Code  7701 
ATTN:  Code  5410,  John  Davli 

• 

Officer-In-Charge 

Naval  Surface  Weapons  Center 

ATTN:  Code  WA501,  Navy  Nuc.  Prgmi.  Off. 


Director 

Strategic  Systems  Project  Office 
ATTN:  NSP 

DEPARTMENT  OF  THE  AIR  FORCE 

AF  Geophysics  Laboratory,  AFSC 
ATTN:  Charles  Pike 
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DEPARTMENT  OF  DEFENSE  COHTRACTORS  (Continued) 


DEPARTMENT  OF  THE  AIR  FORCE  (Continmd) 

AF  Materials  Laboratory,  APSC 
ATTN:  Library 

AF  Weapons  Laboratory,  AFSC 
ATTN:  SUL 
2 cy  ATTN:  NTS 

2 cy  AHN:  DYC 

Headquarters,  USAF/RD 
ATTN:  RDQSM 

Conmander 

Rom  Air  Developmnt  Center,  AFSC 
ATTN:  Edward  A.  Burke 

SAMSO/OY 

ATTN:  DYS 

SAMSO/MN 

ATTN:  MNNH 
ATTN:  MNNG 

SAMS0/5K 

ATTN:  SKF 

SAMSO/XR 

ATTN:  XRS 

Conmander  In  Chief 
Strategic  Air  Comnand 
ATTN:  XPFS 

ATTN:  NRI-STINFO,  Library 

DEPARTMENT  OF  ENERGY 

University  of  California 
Lawrence  Livermore  Laboratory 

ATTN:  Tech.  Info.  Dept.,  L-3 

Los  Alamos  Scientific  Laboratory 

ATTN:  Doc.  Control  for  Reports  Library 

Sandia  Laboratories 
Livermore  Laboratory 

ATTN:  Doc.  Control  for  Theodore  A.  Dell  In 
Sandia  Laboratories 

ATTN:  Doc.  Control  for  3141,  Sandia  Rpt.  Coll. 
OTHER  GOVERNMENT  AGENCY 
NASA 

Lewis  Research  Canter 

ATTN:  N.  J.  Stevens 
ATTN:  Library 
ATTN:  Carolyn  Purvis 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corporation 

ATTN:  Frank  Hal 
ATTN:  Julian  Relnhelmer 
ATTN:  V.  Josephson 
ATTN:  Library 


AVCO  Research  & Systems  Group 

ATTN:  Research  Lib.,  A830,  Rm.  7201 

The  Boeing  Company 

ATTN:  Preston  Qeren 

University  of  Cal Ifornia  at  San  Diego 
ATTN:  Sherman  Da  Forest 

Computer  Sciences  Corporation 
AHN:  Alvin  T.  Schiff 

Dr.  Eugene  P.  DaPlomb 

ATTN:  Eugene  P.  DePlomb 

Dikewood  Industries,  Inc. 

ATTN:  Tech.  Library 
AHN:  K.  Lee 

EQiiG,  Inc. 

Albuquerque  Division 

ATTN;  Technical  Library 

Ford  Aerospace  A Conmunicatlona  Corp. 

ATTN:  Donald  R.  McMorrow,  MS  Q30 
ATTN:  Library 

General  Electric  Company 

Space  Division 

Valley  Forge  Space  Center 

ATTN:  Joseph  C<  Peden,  VFSC,  Rm  4230M 

General  Electric  Company 
TEMPO'Center  for  Advanced  Studies 
ATTN:  William  McNamara 
ATTN;  DASIAC 

Hughes  Aircraft  Cooipany 

ATTN:  Tech.  Library 

Hughes  Aircraft  Company,  El  Ssgundo  Site 
ATTN:  Edward  C,  Smith,  MS  A620 
ATTN:  William  U.  Scott,  MS  AlOBO 

Institute  for  Defense  Analyses 
ATTN:  IDA,  Librarian 

IRT  Corporation 

ATTN:  Dennis  Swift 
ATTN:  Technical  Library 

JAYCOR 

ATTN:  Library 
ATTN:  Eric  P.  Wenaas 

JAYCOR 

ATTN:  Robert  Sullivan 

Johns  Hopkins  University 
Applied  Physics  Laboratory 

ATTN:  Peter  E.  Partridge 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Kaman  Sclancai  Corporation 

I ATTN:  U.  Foiter  Rich 

ATTN:  Jarry  I.  Lubtll 
ATTN:  Library 

LockhMd  Nitallti  A Spaca  Co.,  Inc. 

ATTN:  Dipt.  85-85 

NcDonntll  Douglas  Corporation 

ATTN:  Stanity  Schnaldar 

Mission  Rastarch  Corporation 

ATTN:  Conrad  L.  Longmlra 
ATTN:  Rogar  Stattnar 
ATTN:  N.  J.  Carron 
6 cy  ATTN:  Tech.  Library 

Mission  Rasaarch  Corporat1on-Sa-i  Diago 
AHN:  Library 

I ATTN:  Victor  A.  J.  Van  Lint 

I Rockwall  Intarnatlonal  Corporation 

AHN:  Tachnical  Library 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

RAO  Associates 

ATTN:  Technical  Library 
ATTN:  Leonard  Schlasslngar 

Science  Applications,  Incorporated 
ATTN:  Will  Iasi  L.  Chadsay 

Spire  Corporation 

ATTN:  Rogar  G.  Little 

SRI  International 

ATTN:  Library 

Systans,  Science  and  Softwara,  Inc. 

ATTN:  Technical  Library 
ATTN:  Andrew  R.  Wilson 

TRW  Dafansa  A Spaca  Sys.  Group 

ATTN:  Tech.  Info.  Cantar/S-1930 
?.  cy  AHN:  Robert  M.  Webb,  Rl-2410 
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